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Abstract 

We derive and employ a new semiclassical model to interpret the mechanism of 
the current-induced stress. We assume the metal clusters in spherical jellium model as 
hard sphere and use spherical flow model in classical fluid mechanics and quantum 
viscosity to calculate the cluster radius. The results indicate that there might be 
large-group transportation but not only single-particle transport in electromigration. 
Two new kinds of current-induced stress are derived based on present model, which 
are quantum density stress and quantum viscous stress. These two kinds of quantum 
stress provide new ways to estimate the current- induced stress. Our model is also 
suitable for the motion of nanoisland and void. Our model provides new sight into 
applying hydrodynamics to electromigration. 



Introduction. — Passage of large current in a metal can generate preferential 
diffusion of atoms in one direction. This preferential diffusion of atoms builds up 
compressive stresses near the anode due to the addition of atoms, whereas tensile 
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stresses are found near the cathode due to the removal of material. " The mechanical 
stress leads to grain boundary and interface diffusion, which is known to be driven by 
electron wind force. 4 " 6 The effect of current-induced stress in atomic and molecular 

7 17 

wires has been studied in recent studies. " Previous experiments have found that 
electromigration can have a profound influence upon the development of surface 
morphology and material plasticity. " The mechanism of electron wind force had 
also been broadly studied. " However, some important properties of conducting 
electrons failed to be considered, such as the quantum viscosity. Also, some 
phenomenon cannot be well explained, such as the grain boundary drag effect, the 
motion of nanoislands. 17 

In this article, we provide a new way to interpret the mechanism of 
current-induced stress. We assume the conduction electrons as electron fluid. We 
mainly use the spherical jellium model, ' the quantum Bernoulli equation , the 

71 77 

quantum viscosity , and the classical spherical flow model in fluid mechanics . 
According to the spherical jellium model, nearspherical metal clusters with specific 
number of valence electrons, e.g., 8, 18, 20, 34, 40, 58 are very strong because of their 
enclosed electronic shells. ' These clusters are called magic-number clusters due to 
their exceptional stability. ' The jellium model had been used to analyze the grain 
boundary drag effect. " The grain boundary is regarded as a slab like domain, which 
presents a repulsive potential to the electrons as a shield. Different from the quantum 
mechanical method KKR used in reference 24-26, we assume there exists spherical 
metal clusters in electromigration and assume them as hard sphere. We firstly use 
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quantum Bernoulli equation combined with spherical flow model to analyze the 
current-indcued stress caused by electron density gradient, then we analyze the 
current-indcued stress caused by quantum viscosity The calculated results agree well 
with measured values. Then, we calculated the cluster radius and the number of 
electron per cluster and the calculated results corresponds well with magic numbers, 
which indicates that there might be large group of atoms transportation during 
electromigration but not just single-particle-transport. 

Quantum density stress (QDS). — Electron fluid is the relative reference frame, 
and it moves with velocity V {V = v, -v at0 ) , v ato is the spherical cluster velocity 

which is in the current direction, v e is the free electron drift velocity. On the sphere 

3 
surface, the distribution of relative velocity is v r=a = — V smO , where 9 is the angle 

9 1 

between the radius and x-coordinate. Average the horizontal component of velocity 
of the left side of sphere, we arrive at 

v=- l-Vsm 2 0d0 = -V = -v (1) 

i 

3 
Similarly, the velocity of the other side can be got as — v, . Assuming that A is at 

the left side of dislocation and B is at infinity. The mean horizontal velocity of the 

electron gas at A is v A , the horizontal velocity of the electron gas at B is v e , which is 

the electron gas drift velocity. The pressure of the electron gas at A is P A . The 

pressure of the electron gas at B is P B = P and p B = p , where P is the mean 

electron fluid pressure, p is the mean electron fluid mass density. Each side of the 

dislocation atom is assumed to meet Quantum Bernoulli equation with different 



Generated by Foxit PDF Creator© Foxit Software 
http://www.foxitsoftware.com For evaluation only. 



density, respectively. We use the Bernoulli equation for electron fluid developed by M. 



Di. Ventra 



v 2 P V. 



+ — + -^- = constant (2) 

2 P e m e 

where p e = p I m e , p e is the electron distribution density, m e is the electron mass. 
P = pE f . For conduction electrons, only those very close to the Fermi energy can 
participate in the conduction process. Therefore, their pressure is P = pE f , where 
E , is Fermi level. The pressure equation for electron gas at A is obtained as: 

2 2 

P A =^P+p A e - A (3) 



P 



Similarly, for the other side C 



PC D , „ V e ~ V C 



P(= ^ P+Pc iL_^L (4) 

P 2 



3 
where v c = v A = — v e . According to D'Alembert paradox, the pressure and velocity of 

both sides should be the same. However, the electron gradient between two sides of 

the cluster in the dislocations line leads to the quantum density stress. Then, we obtain 

* ( 3 Y 

° QDS =P A -Pc =^P + Ap e y CJ =^P + ^-Apmv e 2 (5) 
P 2 p e 32 

For the second part of equation (5), it's in fact the kinetic energy transfer from 



electrons to ions. The electron fluid flows with velocity v, and transfers E f to ions. 

1 2 11 

Therefore, — my, should be replaced by —E f , where — is the factor on average. 

2 2 2 

Then equation (5) can be rewritten as 

° EDF = ^pEf + ^ ^pE f = 1 .22Ap e E f (6) 

Ap e is the electron distribution density gradient between two sides of clusters. 

4 
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The distribution of electron is significantly affected by various factors such as charged 
vacancies and congregated solute atoms. "* Between the two sides of dislocation, 
under the effects of charged vacancies, a vivid contrast of electron distribution density 
is formed. The electron density gradient has been observed in previous researches. 
The constrictions in current flow path can increase local density up to about a factor 
of 2. Also, atomic relaxation at step edges redistributed local charge density, ' 
with enhancement up to a factor of 10 on Cu(532) near kink sites. In reference 36, 
using DFT calculation, the charge density gradient of Ag can be represented as: 

l%> e dx = l(fi L -ti R '%p L (fi tq ,r)-p x (ti tq ,r)]flx/2 (11) 

And the simulation results indicates that could be of 4X10" 27 m" 3 on chain and 
dislocation lines, which is directly suitable for present model. The conduction electron 
density p e for single-crystal Ag is 5.85 x 10 28 m" 3 , Fermi energy is 5.48eV, therefore, 
according to equation (6), cr EDF = 4.25Gpa compared with 2GPa which is the 
current induced stress of Pt thin films. Equation (6) contains two parts. The left part 
derives from the Fermi pressure, while the other part derives from the quantum kinetic 
energy transfer. We use classical fluid mechanics to deduce the basic equations for 
EDF Then, we further transform it into quantum mechanical forms, and the calculated 
results demonstrate the success for the transformation we have made. 

Quantum viscous Stress (QVS). — In the first part, we consider the electron fluid 
as ideal fluid and analyze the effect of electron density difference/gradient. In this part, 
we assume the electron fluid as viscous fluid and analyze the viscosity-induced 
pressure. The d.c. viscosity of homogeneous electron fluid could be obtained by AK's 
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approach 



20 
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1.81x10 
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n(l + 2f s ) n 

4 



8^ 



r +r 



(7) 



where r s -ar s ln and a = (4/ 9k) =0.52. For analyzing the effect of quantum 
viscosity, we mainly use hydrodynamics. In classical spherical flow model, under 
stokes linearization approximation, the 2-dimensional pressure distribution on the 



sphere surface could be represented as 



20 



P = R _m^ = PeE 3V./IX 



2R l 



2R 1 



(8) 



where R is the radius of the sphere. The mean pressure on the left side is 

P = p E , . h — . Similarly, the mean pressure of the right side is P = p E, — . 

/ r» f 4R v 6 , ye / 4R 

The quantum viscous stress is thus the difference between the two sides. 



3v eJ u 3jh 



a QVS .-, n .-, .„...(> 
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(9) 



The quantum viscous stress, in fact, could be influenced by some other carrier 
scattering mechanisms or lattice vibration. Therefore, we should add a constant K into 
equation (9) to represent the effective quantum viscous stress. When K=l/2, we find a 
very interesting result. 

TABLE 1 . j is the current density, T is the temperature, a is the previous results 
of current-induced stress. N is the number of electrons inside each characteristic 



sphere and N 



4tt(R^ 



v * J 



Z where Z is the number of free electrons per atom. / 



is the average distance between atoms, R is the radius of metal cluster. 
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System 
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r 
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/ 
10 m 


R 

(K=l/2) 
10- 10 m 


N 


Al conductor line 
(measured) 


0.4 


534 


210 


0.34 


2.55 


3.64 


36.55 


Al conductor line 
(calculated) 


1 


554 


510 


0.34 


2.55 


3.49 


32.21 


Pt thin film (measured) 


2.5 


350 


2000 


0.38 


2.47 


3.25 


19.08 


Sn solder joint 
(calculated) 


0.14 


360 


155 


0.37 


3.00 


2.52 


9.93 



In table 1, the number of electrons inside each cluster of Al, Pt and Sn 
corresponds to the magic number 34, 20, 8. More surprisingly, the calculation 
indicates that the 3 different kinds of material, at different temperature, under different 
current density, in 4 independent experiments, 5 ' 7 ' 9 ' 10 have the same order of magnitude 
of the cluster radius. This could be a strong corroboration to our assumption which is 
based on spherical jellium model. According to spherical jellium model, nearspherical 
metal clusters with magic number are exceptionally strong because of their closed 

1 O 1 Q 11 

electronic shells. ' ' This might indicate that the electromigration transport should 
not only be single-particle transport, a cluster of atoms might be transported altogether. 
It is worth noting that Stokes approximation is satisfied when the Reynolds number 
R e = v e L I jupm e is much smaller than 1 . Therefore, we should calculate the value 

of Reynolds number of electron fluid. According to the value in Table 1, the Reynolds 
number is much less than 1, even when the characteristic size L of the system is 

10 7 m. Besides, equation (9) could only provide a very rough estimation of cluster 
radius, which mainly determines the value of L , a more accurate and well 



correlated approach is needed to get the accurate value of cluster radius. 
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The electron viscous fluid could also be applied to explain the motion of 
nanoislands. Chenggang Tao et al observed the current-biased displacement of 

1 7 

monatomic islands of several nanometers on single-crystal Ag. It was found that the 
velocity of the islands varies inversely with radius, which is contradictory to the 
prediction of pure attachment/detachment mechanism. However, this phenomenon 
can be well explained by our model. Assuming the nanoisland is a disk, then the drag 
force of the disk could be represented as F disk = KVR Q , 21 where K is a constant in this 
system, the island area density is p js , the island radius is R Q , the island drift velocity 



is v . 



Then, we can get 



d^ =K ^A (10) 



K 

v e -Ae p '^ =v (ii) 
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K 
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Because v =0 when / = , then A = v e , therefore 



(12) 

(13) 

According to equation (13), the island drift velocity v varies inversely with 
radius R Q , which well suits the experimental results. Additionally, according to 
equation (13), the island drift velocity will never be larger than the electron drift 
velocity. This phenomenon is thus explained in a completely different way from the 

i n 

diffusion mechanism. An important observation was made by K.C. Chen et al. They 
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found that the grain boundary slowed down electro migration rate by one order of 
magnitude. According to (6), the electromigration direction might be opposite to the 
electron density gradient in grain boundary, which leads to the drag effect. 

In reference 37,38, a solution has been obtained for the motion of void driven by 
electric potential and gradient stress field. If we consider the contribution of QDS and 
QVS, this chemical potential of the atoms on the interface should be rewritten as 

/u = /u -Q.yK-Z*eU-Q.a'-Q.(7 n +Q.U s (14) 

where ju Q is the reference value of the potential, y s is the interface energy, k 

is the curvature of the interface, Q is the atomic volume, U is the applied electric 
potential, U s is the strain energy density of matrix material, cr n is the normal stress 

acting on the interface. <x' represents the quantum density pressure and quantum 
viscous pressure. The atomic flux along the interface is given by 
DAdfi da 



J= — 



die , da da n dU 

y. Ze£;t + — - (15) 

I J S\ (AAA V / 

OS OS OS OS 



ClkT ds KT 

where D j is the interface diffusivity, 5 t is the thickness of the interface, kT 
has the usual meaning, s is the arc length of the interface, E t = -V S U . The mass 
conservation requires that dJ/dy = V/Q. , then integrate it from y to 
y = dJ{\ - m) / (l + m) , where m is the shape factor, we get 

+ C (16) 



Vl + m 



where the constant C can be determined by boundary condition, m is the 
shape factor, then 



c- DS - 



I I 



kT 



V^+V.a' + ^-^l (17) 

ds ds J 
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Symmetry requires that the flux at y = must be zero, thus 

( 



CID& \\ + m 



_»_ ^ , dcr n dU 

eZ E : + V„cr - + 

ds ds 



i s 



(18) 



dkT \ 1 — m 

According to equation (18), quantum density stress contributes to the motion of 
voids, cr' helps explain that the variation of the void velocity should result from the 
variation of electron density distribution. ' Also, QVS of cr' in equation (15) helps 
explain the rate of grain boundary diffusion increasing with the current density 7 

Conclusion. — The current-induced force had been called "electron wind force" 
for a long time, however, it was merely a concept. In present model, this force is 
finally demonstrated as a real electron wind. Our model assumes the electron gas as 
Newton fluid, which is a very bold assumption and completely different from 
previous researches. However, the calculated results and well explanation of relative 
phenomenon strongly corroborates our assumption. We calculate the electron force 
using quantum viscosity, Fermi gas, Stokes drag pressure and Quantum Bernoulli 
formula, and the calculated value agrees well with experimental value. Our model is 
also suitable for the motion and force analysis of materials of various scales, such as 
atoms, nanoislands, dislocations and voids. 
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